The ground-state spectroscopic constants of Be2 revisited 
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^ • Abstract 

g: 

Q , Extensive ab initio calibration calculations combined with extrapolations to- 

^, 
^ I wards the infinite-basis limit lead to a ground-state dissociation energy of 

O 



X 



Be2, De=944ib25 cm~^, substantially higher than the accepted experimental 
value, and confirming recent theoretical findings. Our best computed spectro- 
scopic observables (expt. values in parameters) are G(l) — G(0)=223.7 (223.8), 



^ ; G(2) - G(l)=173.8 (169±3), G(3) - G(2)=125.4 (122±3), and So=0.6086 

0^ 

(0.609) cm ; revised spectroscopic constants are proposed. Multireference 
O '. 

^^ [ calculations based on a full valence CAS (4/8) reference space suffer from 

o\ , 

Q^^ , an unbalanced description of angular correlation; for the utmost accuracy, 

O I the (3s, 3p) orbitals should be added to the reference space. The quality of 

t>^| computed coupled cluster results depends crucially on the description of con- 

*^' nected triple excitations; the CC5SD(T) method yields unusually good results 

because of an error compensation. 

I. INTRODUCTION 

Despite the small size of the beryllium dimer, Be2, a correct computational description 
of its X ^S+ ground state has long been considered as one of the most challenging problems 
in quantum chemistry. |1|] Intuitively one would expect a purely repulsive potential between 
two closed-shell singlet atoms — or perhaps a shallow van der Waals-like minimum — and 



in fact the Hartree-Fock potential is purely repulsive. However, the small (2s) — (2p) gap 
in atomic beryllium complicates the picture, and when angular correlation is admitted, a 
tightly bound molecule is in fact found due to an avoided crossing between (2s)^ + (2s)^ and 
(2s)^(2p^)^ + {2sy{2pzy curves. As a result, the wave function is strongly biconfigurational, 
and in fact an active space of at least four orbitals (the abovementioned plus (25)^(2^2)^ + 
(2s)^ and (2s)^ + (2s)^(2p^)^) is required to obtain a qualitatively correct potential curve [H. 
The Hartree-Fock limit potential is purely repulsive, and early coupled cluster with all 
double excitations (CCD) calculations |^ found only a shallow van der Waals-like minimum. 
Multireference configuration interaction studies 0,§] on the other hand predicted a tightly 
bound minimum, as did (with a highly exaggerated binding energy) a pioneering density 
functional study [|]. These conclusions were corroborated in 1983 by a valence FCI (full 
configuration interaction) study 0, and in the next year, Bondybey and English p[ re- 
ported the first experimental observation. Bondybey subsequently reported i?e=2.45 A 
and the first four vibrational quanta 223.2, 169.7, 122.5, and 79 cm~^; assuming a Morse 
potential, he suggested a dissociation energy of 790±30 cm^^. Petersson and Shirley (PS) 
]TI| , following ab initio calculations of their own, re-analyzed the experimental data in terms 
of a Morse-|-l/-R^ potential and suggested an upward revision to -De=839±10 cm^^. Recent 
high-level calculations suggest even higher binding energies: for instance, Starck and Meyer 
||TT| (SM), using MRCI (multireference configuration interaction) and a core polarization 
potential (CPP) found L'e=893 cm^^ as well as re=2.4485 A, while MR-AQCC (multirefer- 
ence averaged quadratic coupled cluster |jl2|]) calculations by Fiisti-Molnar and Szalay ]13[ 
(FS) established L'e=864 cm~^ as a lower bound. R0eggen and Almlof (RA) [|l|] carried out 
extensive calibration calculations with an extended geminal model and gave 841±18 cm~^ 
as their best estimated binding energy. Evangelisti et al. (EBG) [0] carried out valence- 
only FCI calculations in a [6s5p3d2flg] basis set, and concluded that inner-shell correlation 
must contribute substantially to the binding energy since their value (an exact valence-only 
solution within this large basis set) was still appreciably removed from experiment. This 
conclusion was confirmed by an all-electron FCI in a small [9s2pl(i] basis set (which still 



involved in excess of 10^ determinants) 

Part of the uncertainty in the best theoretical values resides in the fact that the basis 
sets used, while quite large, are still finite. Convergence of angular correlation is known 
to be excruciatingly slow, with an asymptotic expansion in terms of the maximum angular 
momentum / that starts at /~^ for contributions of individual angular momenta and at l^^ for 
overall /-truncation error [|l6l. Recently /-extrapolations have been proposed |17| , [18| which 



permitted the calculation of total atomization energies of small polyatomic molecules with 
mean absolute errors as low as 0.12 kcal/mol. Among other applications, this method made 
possible a definitive re-evaluation [|l^ of the heat of vaporization of boron from a calibration 
quality calculation on BF3. 

In the present work, we apply this method to the dissociation energy of Be2. It will be 
shown that the valence-only basis set limit is in fact as large as 875±10 cm^^, and the overall 
De as large as 945±20 cm~^. 

II. METHODS 

The multireference and FCI calculations, as well as those using the CCSD(T) |^0[ coupled 
cluster method, were carried out using a prerelease version of MOLPR097 Q running on an 
SGI Origin 2000 minisupercomputer at the Weizmann Institute of Science. Calculations with 
other coupled cluster methods were carried out using ACES 110 running on a DEC Alpha 



^ MOLPRO 97.3 is a package of ab initio programs written by H.-J. Werner, and P. J. Knowles, 
with contributions from J. Almlof, R. D. Amos, A. Berning, D. L. Cooper, M. J. O. Deegan, 
A. J. Dobbyn, F. Eckert, S. T. Elbert, C. Hampel, R. Lindh, A. W. Lloyd, W. Meyer, A. 
Nicklass, K. A. Peterson, R. M. Pitzer, A. J. Stone, P. R. Taylor, M. E. Mura, P. Pulay, M. 
Schiitz, H. Stoll, and T. Thorsteinsson, 

2 J. F. Stanton, J. Gauss, J. D. Watts, W. Lauderdale, and R. J. Bartlett, (1996) ACES II, an 
ab initio program system, incorporating the MOLECULE vectorized molecular integral program by 



workstation. 

Most basis sets used belong to the correlation consistent polarized valence n-tuple zeta 
(cc-pVnZ) family of Dunning [^. The cc-pVDZ, cc-pVTZ, cc-pVQZ and cc-pV5Z basis 
sets are [3s2pl(i], [4s3p2(il/], [5s4:p3d2flg], and [6s5p4:d3f2glh] contractions, respectively, 
of (9s4pl(i), {Ils5p2sld) , {12s6p3d2flg), and {14:s8p4:d3f2glh) primitive sets. For assessing 
inner-shell correlation effects, we used the core correlation basis set of Martin and Taylor 
P^ : MTvtz and MTvqz denote completely uncontracted cc-pVTZ and cc-pVQZ basis sets. 



respectively, augmented with one tight p, three tight d, and two tight / functions with 
exponents derived by successively multiplying the highest exponent already in the basis set 
with a factor of three. The MTv5z basis set is obtained similarly, but in addition has a single 
tight g function as well. 

III. RESULTS AND DISCUSSION 
A. Valence electron contribution 

For the cc-pVDZ, cc-pVTZ, and cc-pVQZ basis sets, valence-only FCI calculations could 
be carried out. The results at the reference geometry R = 2.45 A are given in Table 1. 

By comparison with CCD, CCSD |^, and CCSDT |2^ results in the same basis sets 
(CCSDTQ being equivalent to FCI for this case), we can partition the valence binding energy 
into contributions from connected single, double, triple, and quadruple excitations as well 
as investigate their basis set convergence. As previously noted by Sosa et al. [^ in small 
basis sets, no covalent binding is seen at the CCSD level; they found CCSDT-l{a,b} and 
CCSDT-2 to display only a shallow ripple, while CCSDT-4 slightly exaggerates the potential 
well and full CCSDT is slightly above the FCI result. These conclusions are confirmed here; 
moreover, as the basis set is increased, the CCSDT results closely track the FCI ones, which 



J. Almlof, J. and P. R. Taylor, and a modified version of the ABACUS integral derivative package 
by T. Helgaker, H. J. Aa. Jensen, P. J0rgensen, J. Olsen, and P. R. Taylor. 
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in this case implies that the contribution of connected quadruples to the binding converges 
very rapidly to an estimated basis set limit of 85 cm~^. By contrast, the contribution of 
connected triples is actually substantially larger than the atomization energy itself, and is 
apparently not yet converged with the cc-pVQZ basis set. 

Our attempts to carry out a CCSDT/cc-pV5Z calculation with the available computer in- 
frastructure met with failure. CCSD(T) calculations are an obvious alternative, but are seen 
in Table 1 to on the one hand underestimate the importance of connected triple excitations, 
and on the other hand to display considerable basis set dependence in the difference with 
full CCSDT (hence making it a poor candidate for extrapolation). The difference between 
CCSD(T) and CCSDT starts at fifth order in perturbation theory; in the method alterna- 
tively known as CCSD+T(CCSD)* [|6[ and, in Bartlett's recent notation [0, CC5SD(T), 



the missing Er^j-j^ term is included quasiperturbatively at a computational expense scaling as 
'^occ^virf ^s seen in Table 1, CC5SD(T) slightly overestimates the connected triple excita- 
tions contribution but does so in a highly systematic manner, the difference being constant 
between 38 and 40 cm~^. Because of an error compensation with neglect of connected 
quadruple excitations, it is actually the one single-reference method short of full CI that 
we find to be closest to the exact solution. In short, it is the ideal candidate for basis set 
extrapolation. 

The CCSD+TQ(CCSD)* or CC5SD(TQ) method, which includes the leading contribu- 
tion of connected quadruple excitations in a similar fashion, appears to seriously overestimate 
it, and we have not considered it further. 

Basis set superposition error for the valence electrons was considered using the standard 
counterpoise (CP) correction [^. In the present case, it drops from 36 cm~^ (cc-pVDZ) over 



24 (cc-pVTZ) to 6 cm^^ for the cc-pVQZ basis set, and a paltry 3.5 cm^^ for the cc-pV5Z 
basis set. 

From the FCI/cc-pV{D,T,Q}Z results, we may attempt extrapolation, either from the 
uncorrected D^ values (assuming that the extrapolation will absorb BSSE which strictly 
vanishes at the basis set limit) or after subtracting the counterpoise correction in each case. 



With a variable-a 3-parameter correction, this leads to basis set hmits of 841 and 859 cm"^, 
respectively. Using the simple A + B/l^ formula |2^ on just the final two results, we obtain 
values of 863 (raw) and 870 (CP-corrected) cm^^. 

It can rightly be argued that the cc-pVDZ basis set is really too small to be involved 
in this type of extrapolation, and that a cc-pV5Z result is essential for this purpose. This 
requires us to estimate an FCI/cc-pV5Z result from the additivity approximation Method/cc- 
pV5Z+FCI/cc-pVQZ-Method/cc-pVQZ. With Method=CC5SD(T), we obtain De(FCI/cc- 
pV5Z)?^818.2 cm~^; 3-point extrapolation yields 881 cm~^ for the raw, and 872 cm~^ for the 
CP-corrected, results as the basis set limit. Using the simple A + B/l^ formula, we obtain 
the alternative results 857 and 873 cm~^, respectively. The fact that the two extrapolations 
yield essentially the same result for the CP-corrected values, as well as that they are in very 
close agreement with the results with the smaller basis sets, is very satisfying. 

It could likewise be argued that in fact the SCF and correlation contributions should be 
handled separately |18|, with an exponential or (/ + 1/2)^^ formula for the SCF contribution 
and an A + B/{1 + 1/2)" oi A + B/l^ formula for the correlation contribution alone. We 
then find that the SCF contribution, with the cc-pV5Z basis set, lies within 3 cm~^ of the 
numerical HF limit; after adding in the basis set limits for the correlation contribution, we 
obtain, after counterpoise correction, 869 cm~^ with the 3-point and 871 cm~^ with the 
2-point formula. 

One further objection would be to the use of even a high-level single-reference method 
for a problem that is intrinsically multireference in character. We have therefore considered 
MRCI (multireference configuration interaction) augmented with the multireference David- 
son correction pO[, MRACPF [^ (multireference averaged coupled pair functional), and 
MRAQCC |T^ (multireference averaged quadruples coupled cluster) methods with a variety 
of active spaces. A 4/4 active space appears to be unsatisfactory for our purposes; hence 
we have considered full- valence CAS(4/8)-ACPF (averaged coupled pair functional ||31[]) and 
CAS(4/8)-AQCC as alternatives. Except for the cc-pVDZ basis set, both methods seem to 
track the FCI results quite closely, with CAS(4/8)-ACPF accidentally coinciding with the 
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FCI results. Again applying the same additivity approximation as above, we obtain esti- 
mated FCI/cc-pV5Z results from these calculations of 821.5 and 819.6 cm~^, especially the 
latter quite close to the CC5SD(T) derived value. 

Interestingly, the CAS(4/8)-ACPF wave function contains a fairly large number of exter- 
nal excitations with fairly high amplitudes, most of them involving excitation into (3p)-type 
Rydberg orbitals. Inspection of the atomic wave function for Be atom revealed that excita- 
tions into the fairly low-lying (3p) orbitals have amplitudes as large as 0.09 (for each of three 
symmetry-equivalent components); since in addition the (3s) orbital is below the (3p) or- 
bital in energy and there appears to be no clear separation between (3s)- and (3p2)-derived 
a orbitals, this suggests a (4/16) active space which spans all molecular orbitals derived 
from atomic (2s, 2p, 3s, 3p) orbitals. External excitations now carry so little weight in the 
wave function that CAS(4/16)-MRCI+Dav, CAS(4/16)-ACPF and CAS(4/16)-AQCC yield 
essentially identical results. Arbitrarily selecting the CAS(4/16)-ACPF result for extrapola- 
tion, we obtain a best estimate of 821.5 cm~^ for the FCI/cc-pV5Z Df.. After counterpoise 
correction, the CAS(4/16)-ACPF derived value leads to a basis set limit value of 885.6 cm~^ 
with the 3-point and 861.4 cm~^ with the 2-point formula. Taking the average of the latter 
two values and the CC5SD(T) derived ones, we finally propose 872±15 cm^^ as our best 
estimate for the valence-only D^- 

As a final remark, let it be noted that the extrapolations in all cases bridge an area of no 
more than 50-70 cm~^; by substituting / = 6 in the extrapolation fomulas, we can estimate 
that calculations with the next large basis set, cc-pV6Z (i.e. [7s6p5d4f3g2hli]), would only 
recover about 20-25 cm^^ of that total. 

B. Inner-shell contribution 

By taking the difference between their computed MRCI results with and without the 
core polarization potential, SM found that inner-shell correlation would add 0.38 in.Eh, or 83 
cm^^, to the atomization energy. RA computed a contribution of (Is) correlation (almost 



exclusively core-valence correlation) of 0.40654 mEh, or 89.2 cm"^. 

Our results for the effect of inner-shell correlation are collected in Table 2. Using the 
MTvtz, MTvqz, and MTv5z basis sets in succession at the CAS(4/16)-ACPF level, we find 
contributions of inner-shell correlation to the binding energy of 82.1, 80.6, and 77.8 cm"^. 
ESSE contributions to the core-correlation contribution (taken as the difference between 
all-electron and valence-only BSSEs in the same basis set) are 3.8, 2.9, and 1.5 cm~^, re- 
spectively, such that the counterpoise-corrected values of 78.3, 77.7, and 76.3 cm~^ appear 
to be quite handsomely converged. 

For comparison, the counterpoise-corrected CCSD(T) results are 75.0, 73.1, and 70.9 
cm~^, while a CC5SD(T) /MTvtz calculation yielded 63.3 cm~^ without counterpoise cor- 
rection. CAS(4/4)-ACPF and CAS(4/8)-ACPF calculation actually yielded small negative 
inner-shell correlation contributions which are clearly an artifact of the reference space. 

We also note that the counterpoise-corrected all-electron CAS(4/16)-ACPF/cc-pV5Z De 
of 882.4 kcal/mol is already higher than the FS number, and in fact near the SM value. 
Indeed, since this level of electron correlation appears to systematically underestimate the 
valence binding energy by 15-16 cm^^ compared to FCI (see Table 1), we can establish 900 
cm~^ as a lower limit to De- 

Adding the best inner-shell correlation energy contribution of 76.2 cm~^ to our best 
valence binding energy, we obtain a best estimate for the all-electron binding energy of 
948±20 cm~^, where the increased error bar reflects the added uncertainty in the inner-shell 
contribution. 

The effect of scalar relativistic effects was gauged from Darwin and mass-velocity terms 
obtained from CAS(4/16)-ACPF/MTvqz calculations by perturbation theory |3^. At —4.0 
cm^^, it is essentially negligible. 

Combining our best estimates for valence, inner-shell, and relativistic contributions, we 
finally obtain a best estimate for De(Be2) of 944 ±25 cm~^, which suggests that the PS value 
for Df, may need to be revised upward by as much as 100 cm^^. 



C. Potential curve 

Computed bond distances r^, harmonic frequencies Ue, and the first three anharmonicities 
UeXe, uJeVe, and ujgZe are collected in Table 3. They were obtained by a Dunham analysis 
on eighth-order polynomials fitted to some 25 computed energies at bond distances spaced 
around the putative minimum with distances of 0.02 A. 

While good fits could be obtained to the CCSD(T) and CC5SD(T) results, attempts 
to fit CAS(4/8)-{MRCI,ACPF,AQCC} curves in the same manner met with failure. No 
such problem was encountered with results based on a smaller CAS (4/4) reference wave 
function: investigation of the CASSCF energies revealed that while the CAS (4/4) curve is 
bound, the CAS (4/8) curve is purely repulsive in the region sampled. Further investigation 
revealed that with increasing r, amplitudes for excitations into (3p) derived Rydberg orbitals 
progressive take on pathological dimensions (as large as 0.35): under such circumstances, 
the noisy character of the CAS(4/8)-ACPF potential curves should not come as a surprise. 
As expected, expanding the reference space to CAS(4/16) eliminates the problem, as well 
as restores a bound CASSCF potential curve. Apparently the (2p) and (3p) orbitals are 
close enough in importance that a balanced reference space requires that they either be both 
included or both excluded. 

From comparing CAS(4/16)-ACPF/cc-pVTZ and FCI/cc-pVTZ spectroscopic constants, 
it is obvious that the former treatment is indeed very close to an exact solution and the 
method of choice for 1-particle basis set calibration. CC5SD(T) yields surprisingly good rg 
and uJe values (in fact agreeing more closely with FCI than CCSDT) but strongly overesti- 
mates the anharmonicity of the curve. Performance of CCSD(T) is fairly poor, although the 
quality of the results is still amazing considering the pathological character of the molecule. 

Extension of the basis set to cc-pVQZ has a very significant effect on the spectroscopic 
constants, with rg being shortened by 0.026 A and cjg going up by 16 cm^^. Further extension 
to cc-pV5Z has a much milder effect, and suggests that convergence is being approached for 
the molecular properties. A + B/l^ extrapolation suggests that further basis set extension 



may affect Ve by a further —0.003 A and increase cUe by another +2 cm"^. 

Ideally, we would have liked to present all-electron CAS(4/16)-ACPF/MTv5z curves 
in order to include inner-shell correlation. Since however a single point in such a curve 
took more than a day of CPU time on an SGI Origin 2000, we have not pursued this 
option further, and have instead contented ourselves with considering the difference between 
CCSD(T)/MTv5z curves with and without constraining the (Is)-like orbitals to be doubly 
occupied. Our results suggest that inner-shell correlation reduces r^ by 0.03 A and increases 
ujf. by 14 cm~^. The spectroscopic constants given as 'best estimate' are obtained by adding 
these contributions to the extrapolated CAS(4/16)-ACPF/cc-pVooZ results, as well as the 
small difference between FCl/cc-pVTZ and CAS(4/16)-ACPF/cc-pVTZ. 

Obviously, given the highly anharmonic nature of the potential surface, a Dunham-type 
perturbation theory analysis is not appropriate. Like in our recent calibration study on 
the first-row diatomic hydrides, we have transformed our 8th-order Dunham expansion and 
computed dissociation energy to a variable-beta Morse (VBM) potential [Q 



Vc = De(l- exp[-z(l + biZ + b2Z^ + ... + &6^^)])' (1) 

in which z = (3{r — rf.)/rf. and the parameters 6„ and (3 are obtained by derivative matching 



as discussed in detail in Ref. |3J]. The one-dimensional Schrodinger equation was then 
integrated using the algorithm of Balint-Kurti et al. I^Sf, on a grid of 256 points over the 
interval [0.2re,3re]. 

The results for the first four vibrational quanta are given in Table |rV|. We have considered 
three potentials. The first two are the uncorrected FCl/cc-pVTZ and CAS(4/16)-ACPF/cc- 
pV5Z potentials; the third one was obtained by substituting our best estimate D^ and re, 
and adjusting (3 such that the best estimate uj^. is matched. (The 6„ remain unchanged from 
the CAS(4/16)-ACPF/cc-pV5Z values.) What this latter approaches in effect assumes is 
that the shape of the CAS (4/16)- ACPF/cc-pV5Z curve is fundamentally sound. 

As expected, the unadjusted FCI/cc-pVTZ potential seriously underestimates the first 
three vibrational quanta because of the strong dependence of De-, uj^-, and r^ on the basis set 
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and the inclusion of inner-shell correlation. CAS(4/16)-ACPF/cc-pV5Z does so to a much 
lesser extent. Our 'best estimate' potential, however, reproduces the fundamental (the only 
transition known with some precision) essentially exactly, and is in good agreement with 
experiment for the next two quanta. Since the VBM form of the potential does not take 
into account long-distance behavior and the fourth quantum lies at 80% of the dissociation 
energy, it is not surprising that the fourth quantum is seriously overestimated. 

Finally, let us turn to the spectroscopic constants derived from our best potential (Table 
5). Our best Ue is in perfect agreement with SM but substantially lower than the Bondybey 
value. Our best UeXe is substantially smaller than both the Bondybey and SM values: 
however, both of the latter were determined phenomenologically as [^(2) — 2G{1) — G'(0)]/2 
and therefore include contributions from higher-order anharmonicities. If we compute the 
same quantity, we obtain perfect agreement with the SM value. While our rotation- vibration 
coupling constant Ue is in very good agreement with the SM calculations, it is substantially 
larger than the Bondybey value. However, it should be noted that the Be2 potential is 
so anharmonic that the series -B„ = B^ — aein + 1/2) + 7e(n -|- 1/2)^ + 5e{n + l/2)'^ + . . . 
cannot be truncated after the linear term; from our best computed spectroscopic constants, 
we obtain i?o=0.6086 cm^^, in perfect agreement with Bondybey's value of 0.609 cm^^ for 
this observable quantity. In short, we argue that our computed r^ = 2.440 A is more reliable 
than the Bondybey value of 2.45o A. 

As a final note, we point out that this revised reference geometry (re=2.440 A) would not 
have affected our calculation of De materially, since the energy difference between R =2.44 
and R = 2 Ah A with our best potential only amounts to 0.4 cm~^. 

IV. CONCLUSIONS 

From an exhaustive basis set convergence study on the dissociation energy of the ground- 
state Be2, we find that the accepted experimental value needs to be revised upward to a best 
estimate of 944 ±25 cm^^. Individual contributions to this value include a valence-only 
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FCI basis set limit of 872±15 cm~^, an inner-shell contribution of 76±10 cm~^, and rela- 
tivistic corrections as small as —4 cm~^. The performance of single-reference methods for 
this molecule is crucially dependent on their treatment of connected triple excitations; while 
CCSD(T) underestimates binding in this molecule, the CC5SD(T) method performs surpris- 
ingly well at a fraction of the cost of full CCSDT. The contribution of connected quadruple 
excitations is small (80 cm^^) and fairly insensitive to the basis set. Accurate multireference 
calculations require an active space which treats angular (2p,3p) correlation in a balanced 
way; a full- valence CAS (4/8) reference does not satisfy this criterion. For the utmost accu- 
racy, a CAS (4/16) reference including the (3s, 3p) orbitals is required, while for less accurate 
work a CAS (4/4) reference is recommended. Our best computed spectroscopic observables 
(expt. values in parameters) are G(l) - G'(0)=223.7 (223.8), ^(2) - G'(l) = 173.8 (169±3), 
^(3) - G'(2) = 125.4 (122±3), and 5o=0.6086 (0.609) cm^^ Our best computed spectro- 
scopic constants represent substantial revisions from the experimentally derived values; in 
particular, the bond length is 0.01 A shorter than the accepted experimental value. 
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TABLES 
TABLE L Convergence of the valence dissociation energy (cm^^) of Be2 as a function of basis 
set and electron correlation treatment 



FCI 


cc-pVDZ 
23.92 


cc-pVTZ 
630.50 


cc-pVQZ 
764.81 


cc-pV5Z 








Difference with FCI 




Actual 
value 


Estimated 










FCI/ce-pV5Z'^ 



SCF 


-2759.54 


-3277.17 


CAS(4/8)-CI+Davidson 


38.85 


-36.27 


CAS(4/4)-ACPF 


35.66 


-56.63 


CAS(4/4)-AQCC 


22.10 


-84.61 


CAS(4/8)-CI+Davidson 


84.77 


35.53 


CAS(4/8)-ACPF 


60.02 


0.22 


CAS(4/8)-AQCC 


43.41 


-23.76 


CAS(4/16)-CI+Davidson 


48.12 


-14.17 


CAS(4/16)-ACPF 


47.94 


-14.78 


CAS(4/16)-AQCC 


48.26 


-15.25 


CCD 


-978.97 


-1245.24 


CCSD 


-943.86 


-1118.72 


CCSDT 


-73.43 


-83.13 


CCSD(T) 


-261.65 


-257.57 


CC5SD(T) 


-39.75 


-38.15 



3396.48 


-2626.13 


770.36 


-49.38 


769.23 


818.61 


-69.67 


747.16 


816.82 


-98.23 


717.18 


815.41 


36.96 


859.67 


822.71 


-0.50 


821.03 


821.53 


-24.71 


794.89 


819.60 


-14.16 


807.49 


821.65 


-15.05 


806.47 


821.52 


-15.61 


805.80 


821.40 


1278.06 


-474.45 


803.61 


1158.37 


-353.16 


805.21 


-84.47 


— 


— 


-241.37 


587.41 


828.78 


-40.07 


778.09 


818.16 



ESSE" 



36.00 



24.37 



6.10 



3.47 



(a) according to FCI/cc-pV5Z ^ Method/cc-pV5Z 

(b) counterpoise method 



FCI/cc-pVQZ - Method/cc-pVQZ 
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TABLE II. Contribution of inner-shell correlation to the dissociation energy (cm 


1) of Be2 


Method e^ correlated MTvtz MTvqz 


MTv5z 



CCSD(T) 



CC5SD(T) 



CAS(4/4)-ACPF 



CAS(4/8)-ACPF 



CAS(4/16)-ACPF 



ESSE (a) 



all 

valence 

difference 

all 

valence 

difference 

all 

valence 

difference 

all 

valence 

difference 

all 

valence 

difference 

all 

valence 

difference 



507.36 
432.34 

75.03 
641.77 
705.06 

63.29 
580.39 
588.89 

-8.50 
679.29 
682.17 

-2.88 
749.67 
667.56 

82.11 
9.43 
5.63 
3.80 



614.26 

541.15 

73.11 



673.94 

676.80 

-2.86 

773.33 

779.23 

-5.90 

845.37 

764.81 

80.56 

7.22 

4.36 

2.86 



661.56 

590.68 

70.88 



823.22 

811.44 

-11.78 

886.56 

808.71 

77.85 

4.06 

2.52 

1.54 



(a) on CAS(4/16)-ACPF values 
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TABLE III. Convergence of computed bond distance (A) and vibrational spectroscopic con- 
stants (cm~^) of Be2 as a function of basis set and electron correlation treatment 



Method 




Basis 


re 


UJe 


LUeXe 


UJeVe 


IJJeZe 


CC5SD(T) 




cc-pVDZ 


2.5736 


187.4 


33.175 


-4.937 








cc-pVTZ 


2.50f2 


230.8 


23.198 


-1.179 


-0.116 






cc-pVQZ 


2.4745 


245.3 


21.825 


-1.072 








cc-pV5Z 


2.47f8 


247.5 


21.367 


-0.959 




CCSD(T) 




MTvtzALL" 


2.4829 


229.2 


24.689 


-1.603 








MTvtzVAL'' 


2.5145 


214.3 


25.919 


-1.871 








difference 


-0.0316 


14.9 


-1.230 


0.268 








MTvqzALL 


2.4685 


241.5 


22.987 


-1.316 








MTvqzVAL 


2.4986 


227.2 


23.905 


-1.514 








difference 


-0.0301 


14.3 


-0.918 


0.198 








MTv5zALL 


2.4652 


243.5 


22.721 


-1.214 








MTv5zVAL 


2.4950 


229.5 


23.482 


-1.335 








difference 


-0.0298 


14.0 


-0.761 


0.120 




FCI 




cc-pVDZ 


2.5598 


193.9 


31.174 


-4.082 








cc-pVTZ 


2.5021 


234.3 


22.383 


-1.071 


-0.097 


CAS(4/16) 


-ACPF 


cc-pVTZ 


2.5041 


232.4 


22.639 


-1.103 


-0.093 






cc-pVQZ 


2.4781 


246.7 


21.325 


-1.011 


-0.084 






cc-pV5Z 


2.4750 


249.1 


20.856 


-0.905 


-0.061 






cc-pVooZ (c) 


2.4718 


251.7 


20.365 


-0.793 


-0.037 


Best estimate (d) 


2.4397 


267.9 


19.191 


-0.563 


-0.042 




Bondybey 


i 


2.45o 


275.8 


26.0 








SM [[i| 




2.4485 


268.2 


24.9 









(a) all electrons correlated 

(b) only valence electrons correlated 

(c) extrapolated according to A + B/P 

(d) CAS(4/16)-ACPF/cc-pVcx)Z+[FCI/cc-pVTZ-CAS(4/16)-ACPF/cc-pVTZ] 

+ [CCSD(T)/MTv5zALL-CCSD(T)/MTv5zVAL] 
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TABLE IV. Computed and observed vibrational energy level differences (cm ^) for the X ^E"*" 
state of Be2 



ZPE 




110.6 


G(l)- 


-G(0) 


185.4 


G(2)- 


-G(l) 


125.5 


G(3)- 


-G{2) 


72.2 


G(4)- 


-G(3) 


75.3 



118.5 
204.2 
153.5 
109.5 
99.2 



127.9 
223.7 
173.8 
125.6 
106.9 



223.8 
169 
122 
79 



125 

218 
168 
112 
67 



213 

167 
122 

78 



FCI/ CAS(4/16)-ACPF/ best Expt. H |0| lTl| 

cc-pVTZ cc-pV5Z (a) 



124.8 
218.4 
168.6 
112.1 
69.4 



1 



221.0 
162.9 
94.2 

54.7 



(a) from FCI/cc-pVTZ potential in form eq.(l), but with best estimate re, -De, and u^ 

substituted according to /3new//3old = ^e,ne^re,r.e^\[D^s./^e,o\dre,o\d\[D^^ 



TABLE V. Potential function parameters in eq.(l) and mechanical spectroscopic constants of 
Be2 with this potential. All values in cm~^ except /3 and the 6„, which are dimensionless 



best potential 



calculated 



Bondybey 



SM [11 



De=944.0 (a) 
re=2.439685 
/3=5.499750 
fei= 0.019920 
fea =-0.048391 
63 =-0.016734 
64= 0.000693 
65= 0.001938 
b(,= 0.000324 



yoo=-0.788 

(.Je=267.93 

u)eXE=2a.mi (d) 

LJeJ/e=-0.827 

u)^Ze=-0.0b2 

Se=0.62853 

Bo =0.60863 

ae=0.03787 (b) 

7e =-0.00361 

(5e=-0.00050 

De=13.84xl0-"5 (c) 

/3e=3.48xl0-'' 



(a) dissociation energy 

(b) -{Bi - 5o)=0.02904 cm"^ 

(c) quartic centrifugal distortion constant 

(d) [G(2) - 2G(1) - G(0)]/2=24.95 cm'^ 



275.8 
26.0 



0.623 
0.609 
0.028 



14.8 



268.2 
24.9 



0.6213 



0.037 
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